Background

This notebook is for preparing input files for pyclone-vi. The same files with a minor modification will be used as an input for FastClone.

Set up

suppressPackageStartupMessages({
  library(tidyverse)
  library(readxl)
  #library(cDriver) # Calculate CCF, https://github.com/hanasusak/cDriver/
})

Directories and File Inputs/Outputs

# Detect the ".git" folder -- this will be in the project root directory
# Use this as the root directory to ensure proper sourcing of functions 
# no matter where this is called from
root_dir <- rprojroot::find_root(rprojroot::has_dir(".git"))
analysis_dir <- file.path(root_dir, "analyses", "tumor-clone-inference")
input_dir <- file.path(analysis_dir, "input")
data_dir <- file.path(root_dir, "data")
files_dir <- file.path(root_dir, "analyses", "sample-distribution-analysis", "results")
maf_files_dir <- file.path(root_dir, "analyses", "tmb-vaf-longitudinal", "results")
#tmb_input_dir <- file.path(root_dir, "analyses", "tmb-vaf-longitudinal", "input")
#scratch_dir <- file.path(root_dir, "scratch")

# Input files
pbta_file <- file.path(files_dir, "pbta.tsv") # file from add-sample-distribution module
maf_file <- file.path(data_dir, "snv-consensus-plus-hotspots.maf.tsv.gz")
tmb_file <- file.path(maf_files_dir, "tmb_vaf_genomic.tsv")
genomic_paired_file <- file.path(files_dir, "genomic_assays_matched_time_points.tsv") # file from add-sample-distribution module
nautilus_dec_file <- file.path(input_dir, "deceased_samples.xlsx") 
palette_file <- file.path(root_dir, "figures", "palettes", "tumor_descriptor_color_palette.tsv")
#tmb_all_file <- file.path(tmb_input_dir, "snv-mutation-tmb-coding.tsv")


# EXAMPLE 
# But to replace with dir with cns data inferred by CNVkit
# PT_Z4BF2NSB_dir <- file.path(input_dir, "cnvkit_data_example/PT_Z4BF2NSB") 
cns_dir <- file.path(input_dir, "cnvkit_data") 

# File path to results directory
results_dir <-
  file.path(analysis_dir, "results")
if (!dir.exists(results_dir)) {
  dir.create(results_dir)
}

# File path to input directory
pyclonevi_input_dir <-
  file.path(analysis_dir, "results", "pyclone-vi-input")
if (!dir.exists(pyclonevi_input_dir)) {
  dir.create(pyclonevi_input_dir)
}

# File path to input directory
fastclone_input_dir <-
  file.path(analysis_dir, "results", "fastclone-input")
if (!dir.exists(pyclonevi_input_dir)) {
  dir.create(pyclonevi_input_dir)
}


# File path to plot directory
pyclone_plots_dir <-
  file.path(analysis_dir, "plots", "pyclone-vi")
if (!dir.exists(pyclone_plots_dir)) {
  dir.create(pyclone_plots_dir)
}


source(paste0(root_dir, "/figures/scripts/theme.R"))

Load and process data

# Let's count #specimens per sample 
# We will remove any samples with less than 2 specimens as phylogenies require at least 3 taxa
kids_specimens_n_df <- genomic_paired_df %>% 
  select(Kids_First_Participant_ID, Kids_First_Biospecimen_ID, tumor_descriptor) %>% 
  unique() %>% 
  dplyr::count(Kids_First_Participant_ID) %>% 
  dplyr::mutate(kids_specimens_n = glue::glue("{Kids_First_Participant_ID}  (N={n})")) %>%
  dplyr::rename(kids_specimens_number = n) %>% 
  filter(!kids_specimens_number <= 2) %>%

  left_join(genomic_paired_df, by = c("Kids_First_Participant_ID")) %>%
  #left_join(maf_df) %>% 
  filter(!is.na(Chromosome)) 


# Let's confirm and add the number of timepoints per sample
# In case we lost any from filtering hypermutants and high reads/alteration
timepoints_n_df <- kids_specimens_n_df %>% 
  select(Kids_First_Participant_ID, tumor_descriptor) %>% 
  unique() %>% 
  dplyr::count(Kids_First_Participant_ID) %>% 
  dplyr::mutate(kids_timepoints_n = glue::glue("{Kids_First_Participant_ID}  (N={n})")) %>%
  dplyr::rename(kids_timepoints_number = n) %>% 
  filter(!kids_timepoints_number <= 1)

df <- timepoints_n_df %>% 
  left_join(kids_specimens_n_df, by = c("Kids_First_Participant_ID")) %>% 
  write_tsv(file.path(results_dir, "samples_eligible_for_phylogeny.tsv"))

# Number of samples in the `genomic_paired_df` with > 3 biospecimens and >2 timepoints
samples_with_3bs <- print(length(unique(df$Kids_First_Participant_ID)))
[1] 34
samples_with_3bs <- print(unique(df$Kids_First_Participant_ID))
 [1] "PT_04V47WFC" "PT_0DWRY9ZX" "PT_19GCSK2S" "PT_1H2REHT2" "PT_23NZGSRJ" "PT_3KM9W8S8" "PT_75HRTX4S" "PT_8GN3TQRM" "PT_962TCBVR" "PT_99S5BPE3" "PT_9S6WMQ92"
[12] "PT_AQWDQW27" "PT_C1RDBCVM" "PT_CWXSP19D" "PT_D6AJHDST" "PT_DVXE38EX" "PT_EQX0VT4F" "PT_GTHZF21E" "PT_HFQNKP5X" "PT_HHG37M6W" "PT_HJMP6PH2" "PT_K8ZV7APT"
[23] "PT_KBFM551M" "PT_KTRJ8TFY" "PT_KZ56XHJT" "PT_MDWPRDBT" "PT_MNSEJCDM" "PT_N8W26H19" "PT_NJQ26FHN" "PT_NK8A49X5" "PT_S4YNE17X" "PT_TP6GS00H" "PT_Z4BF2NSB"
[34] "PT_ZZRBX5JT"
  
# List with samples eligible for phylogeny
# I added the information about to use for phylogenetic inferences
list_df <- df %>% 
  select(Kids_First_Participant_ID) %>% 
  unique() %>% 
  mutate(somatic_germline_phylogeny = case_when(grepl("PT_KZ56XHJT|PT_KTRJ8TFY", Kids_First_Participant_ID) ~ "yes",
                                        TRUE ~ "not_yet")) %>% 
  write_tsv(file.path(results_dir, "samples_eligible_for_phylogeny_list.tsv"))

# Remove large files
rm(genomic_paired_df, kids_specimens_n_df, tmb_df)

Add Nautilus location for Deceased specimens

nautilus_dec_df <- read_excel(nautilus_dec_file) %>% 
  right_join(df, by = c("sample_id", "aliquot_id", "tumor_descriptor")) %>%
  dplyr::mutate(mutation_id = paste(Kids_First_Participant_ID, tumor_descriptor, Kids_First_Biospecimen_ID, Chromosome, Start_Position, Reference_Allele, Tumor_Seq_Allele2, sep = ":")) #%>% 
  #write_tsv(file.path(results_dir, "nautilus_dec.tsv"))

list_df <- nautilus_dec_df %>% 
  select(Kids_First_Participant_ID, Kids_First_Biospecimen_ID, `Note field from Nautilus of initial parent`) %>% 
  unique() %>% 
  write_tsv(file.path(results_dir, "nautilus_dec_list.tsv"))

Other, maybe to delete


# Make list with samples with >2 biospecimens at Deceased timepoint
dec_n_df <- df %>% 
  filter(tumor_descriptor == "Deceased") %>% 
  select(Kids_First_Participant_ID, tumor_descriptor, Kids_First_Biospecimen_ID) %>% 
  unique() %>% 
  dplyr::count(Kids_First_Participant_ID) %>% 
 # dplyr::mutate(kids_dec_n = glue::glue("{Kids_First_Participant_ID}  (N={n})")) %>%
  dplyr::rename(kids_deceased_bs_number = n) %>% 
  filter(!kids_deceased_bs_number <= 1) %>% 
  write_tsv(file.path(results_dir, "kids_dec_multiple_bs_list.tsv"))

###-----------------------------------------------------
# let's look into PT_3KM9W8S8
#PT_3KM9W8S8_df <- nautilus_dec_df %>% 
#  filter(Kids_First_Participant_ID == "PT_3KM9W8S8",
 #        Kids_First_Biospecimen_ID == "BS_2NQXY528") %>% 
#  select(Kids_First_Participant_ID, Kids_First_Biospecimen_ID, `Note field from Nautilus of initial parent`) %>% 
#  unique() %>% 
#  write_tsv(file.path(results_dir, "nautilus_dec_list.tsv"))
  
#bs_id <- print(unique(PT_3KM9W8S8_df$Kids_First_Biospecimen_ID))

Create pyclone input files

We need to generate the input files according to the method’s template. Phylogenetic methods require at least 2 samples per tumor site (multiregional sampling per anatomical site). Here, we will consider kids samples with more than 2 timepoints with one or more biospecimens. We will compare later differences in samples with single vs multiple biospecimens.

# Create pyclone df for all samples
pyclone_all_samples_df <- nautilus_dec_df %>%
  select(Kids_First_Participant_ID, Kids_First_Biospecimen_ID, cg_id, 
         cg_id_kids, mutation_id, sample_id,
         tumor_descriptor, Chromosome, Start_Position, 
         Reference_Allele, Tumor_Seq_Allele2, t_ref_count, t_alt_count,
         normal_cn, tumor_fraction, mutation_count) %>% 
  
  # change names to match input requirements
  dplyr::rename("ref_counts" = "t_ref_count", 
                "alt_counts" = "t_alt_count",
                "tumour_content" = "tumor_fraction") %>%
  select(Kids_First_Participant_ID, Kids_First_Biospecimen_ID, tumor_descriptor, 
         cg_id_kids, mutation_id, sample_id, ref_counts, 
         alt_counts, normal_cn, tumour_content, mutation_count)

samples_pyclone <- print(unique(pyclone_all_samples_df$Kids_First_Participant_ID))
 [1] "PT_04V47WFC" "PT_0DWRY9ZX" "PT_19GCSK2S" "PT_1H2REHT2" "PT_23NZGSRJ" "PT_3KM9W8S8" "PT_75HRTX4S" "PT_8GN3TQRM" "PT_962TCBVR" "PT_99S5BPE3" "PT_9S6WMQ92"
[12] "PT_AQWDQW27" "PT_C1RDBCVM" "PT_CWXSP19D" "PT_D6AJHDST" "PT_DVXE38EX" "PT_EQX0VT4F" "PT_GTHZF21E" "PT_HFQNKP5X" "PT_HHG37M6W" "PT_HJMP6PH2" "PT_K8ZV7APT"
[23] "PT_KBFM551M" "PT_KTRJ8TFY" "PT_KZ56XHJT" "PT_MDWPRDBT" "PT_MNSEJCDM" "PT_N8W26H19" "PT_NJQ26FHN" "PT_NK8A49X5" "PT_S4YNE17X" "PT_TP6GS00H" "PT_Z4BF2NSB"
[34] "PT_ZZRBX5JT"
# Remove large files
# rm(nautilus_dec_df)

# I will test one HGG dataset for now
# PT_Z4BF2NSB

#PT_Z4BF2NSB_DF <- pyclone_all_samples_df %>% 
#  filter(Kids_First_Participant_ID == "PT_Z4BF2NSB")
                   

Add major_cn, minor_cn columns from cns files

samples_pyclone <- print(unique(pyclone_input$Kids_First_Participant_ID))
[1] "PT_KTRJ8TFY" "PT_KZ56XHJT" "PT_Z4BF2NSB"

Save input files

for (i in 1:length(data_dir) ) { 
  
  # Create sample_name
  sample_name <-  unique(as.character(gsub(".call.cns", "", str_split_fixed(data_dir[i], "/", 13)[,11])))
  #sample_name <-  unique(as.character(gsub(".call.cns", "", str_split_fixed(data_dir[i], "/", 10)[,9])))
  sample_name <- sort(sample_name, decreasing = FALSE)
  print(sample_name)
  
for (x in seq_along(sample_name) ) { 
  
  # Save input file for pyclone-vi
  pyclone_fname <- paste0(pyclonevi_input_dir, "/", sample_name[x], ".tsv")
  print(pyclone_fname)
  
  pyclone_input_subset <- pyclone_input %>%
    filter(Kids_First_Participant_ID == sample_name[x]) %>% 
    #select(-c(Kids_First_Participant_ID)) %>% 
  write_tsv(file.path(pyclone_fname))
  
  
  # Save input file for FastClone
  fastclone_fname <- paste0(fastclone_input_dir, "/", sample_name[x], ".tsv")
  print(fastclone_fname)
  
  fastclone_input_subset <- pyclone_input %>%
    dplyr::rename("var_counts" = "alt_counts") %>% 
    filter(Kids_First_Participant_ID == sample_name[x]) %>% 
    #select(-c(Kids_First_Participant_ID)) %>% 
  write_tsv(file.path(fastclone_fname))
  
  }
}
[1] "PT_KTRJ8TFY"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_KTRJ8TFY.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_KTRJ8TFY.tsv"
[1] "PT_KTRJ8TFY"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_KTRJ8TFY.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_KTRJ8TFY.tsv"
[1] "PT_KTRJ8TFY"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_KTRJ8TFY.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_KTRJ8TFY.tsv"
[1] "PT_KTRJ8TFY"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_KTRJ8TFY.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_KTRJ8TFY.tsv"
[1] "PT_KTRJ8TFY"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_KTRJ8TFY.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_KTRJ8TFY.tsv"
[1] "PT_KTRJ8TFY"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_KTRJ8TFY.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_KTRJ8TFY.tsv"
[1] "PT_KTRJ8TFY"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_KTRJ8TFY.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_KTRJ8TFY.tsv"
[1] "PT_KTRJ8TFY"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_KTRJ8TFY.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_KTRJ8TFY.tsv"
[1] "PT_KZ56XHJT"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_KZ56XHJT.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_KZ56XHJT.tsv"
[1] "PT_KZ56XHJT"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_KZ56XHJT.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_KZ56XHJT.tsv"
[1] "PT_KZ56XHJT"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_KZ56XHJT.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_KZ56XHJT.tsv"
[1] "PT_KZ56XHJT"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_KZ56XHJT.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_KZ56XHJT.tsv"
[1] "PT_KZ56XHJT"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_KZ56XHJT.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_KZ56XHJT.tsv"
[1] "PT_KZ56XHJT"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_KZ56XHJT.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_KZ56XHJT.tsv"
[1] "PT_KZ56XHJT"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_KZ56XHJT.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_KZ56XHJT.tsv"
[1] "PT_KZ56XHJT"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_KZ56XHJT.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_KZ56XHJT.tsv"
[1] "PT_KZ56XHJT"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_KZ56XHJT.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_KZ56XHJT.tsv"
[1] "PT_KZ56XHJT"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_KZ56XHJT.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_KZ56XHJT.tsv"
[1] "PT_KZ56XHJT"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_KZ56XHJT.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_KZ56XHJT.tsv"
[1] "PT_KZ56XHJT"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_KZ56XHJT.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_KZ56XHJT.tsv"
[1] "PT_Z4BF2NSB"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_Z4BF2NSB.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_Z4BF2NSB.tsv"
[1] "PT_Z4BF2NSB"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_Z4BF2NSB.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_Z4BF2NSB.tsv"
[1] "PT_Z4BF2NSB"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_Z4BF2NSB.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_Z4BF2NSB.tsv"
[1] "PT_Z4BF2NSB"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_Z4BF2NSB.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_Z4BF2NSB.tsv"
[1] "PT_Z4BF2NSB"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_Z4BF2NSB.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_Z4BF2NSB.tsv"
[1] "PT_Z4BF2NSB"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_Z4BF2NSB.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_Z4BF2NSB.tsv"
[1] "PT_Z4BF2NSB"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_Z4BF2NSB.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_Z4BF2NSB.tsv"
[1] "PT_Z4BF2NSB"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/pyclone-vi-input/PT_Z4BF2NSB.tsv"
[1] "/Users/chronia/CHOP/GitHub/pbta-tumor-evolution/analyses/tumor-clone-inference/results/fastclone-input/PT_Z4BF2NSB.tsv"

Plot depth coverage

# Read color palette
palette_df <- readr::read_tsv(palette_file, guess_max = 100000, show_col_types = FALSE) %>% 
  mutate(tumor_descriptor = color_names)

# Define and order palette
palette <- palette_df$hex_codes
names(palette) <- palette_df$tumor_descriptor

# Define timepoints
timepoints = c("Diagnosis", "Progressive", "Recurrence", "Deceased", "Second Malignancy", "Unavailable")
for (i in 1:length(data_dir) ) { 
  
  # Create sample_name
  sample_name <-  unique(as.character(gsub(".call.cns", "", str_split_fixed(data_dir[i], "/", 13)[,11])))
  sample_name <- sort(sample_name, decreasing = FALSE)
  print(sample_name)
  
for (x in seq_along(sample_name) ) { 
  
  pyclone_input_subset <- pyclone_input %>%
    filter(Kids_First_Participant_ID == sample_name[x]) %>% 
    select(-c(Kids_First_Participant_ID)) #%>% 
   # mutate(tumor_descriptor = factor(tumor_descriptor),
    #     tumor_descriptor = fct_relevel(tumor_descriptor, timepoints)) %>% 
    #arrange(tumor_descriptor, sample_id)
  
  # Make this reproducible
  set.seed(2023)

  # Define label for plots
  Timepoint <- factor(x = pyclone_input_subset$tumor_descriptor, levels = timepoints)
  
  #######################
  # Create bxp ref_counts
    p <- print(ggplot(pyclone_input_subset, aes(sample_id, ref_counts, color = Timepoint)) + 
                 geom_jitter(width = 0.15, size = 0.7, alpha = 0.6) +
                 ggplot2::geom_boxplot(color = "black",
                              size = 0.25,
                              alpha = 0,
                              coef = 0) + # remove whiskers
                 theme_Publication() + 
                 scale_color_manual(values = palette,
                                    breaks = sort(names(palette))) +
                 #rotate() +
                 theme(axis.text.x = element_text(angle = 90)) +
                 stat_summary(fun.y=mean,shape=1,col='black',geom='point') +
                 labs(title = sample_name[x],
                      x = "sample_id",
                      y = "ref_counts",
                      color = "Timepoint"))
    
    # Save the plot
    ggsave(filename = paste0(sample_name[x], "-ref_counts.pdf"), 
         path = pyclone_plots_dir, 
         width = 6, 
         height = 5, 
         device = "pdf", 
        useDingbats = FALSE)
    

  #######################
  # Create bxp alt_counts
    p <- print(ggplot(pyclone_input_subset, aes(sample_id, alt_counts, color = Timepoint)) + 
                 geom_jitter(width = 0.15, size = 0.7, alpha = 0.6) +
                 ggplot2::geom_boxplot(color = "black",
                              size = 0.25,
                              alpha = 0,
                              coef = 0) + # remove whiskers
                 theme_Publication() + 
                 scale_color_manual(values = palette,
                                    breaks = sort(names(palette))) +
                 #rotate() +
                 theme(axis.text.x = element_text(angle = 90)) +
                 stat_summary(fun.y=mean,shape=1,col='black',geom='point') +
                 labs(title = sample_name[x],
                      x = "sample_id",
                      y = "alt_counts",
                      color = "Timepoint"))
    
    # Save the plot
    ggsave(filename = paste0(sample_name[x], "-alt_counts.pdf"), 
         path = pyclone_plots_dir, 
         width = 6, 
         height = 5, 
         device = "pdf", 
        useDingbats = FALSE)
  
  
}
}
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"

Total number of mutations across timepoints and biospecimen sample per Patient case

for (i in 1:length(data_dir) ) { 
  
  # Create sample_name
  sample_name <-  unique(as.character(gsub(".call.cns", "", str_split_fixed(data_dir[i], "/", 13)[,11])))
  sample_name <- sort(sample_name, decreasing = FALSE)
  print(sample_name)
  
for (x in seq_along(sample_name) ) { 
  
  pyclone_input_subset <- pyclone_input %>%
    filter(Kids_First_Participant_ID == sample_name[x]) %>% 
    select(-c(Kids_First_Participant_ID))
  
  # Define parameters for function
  ylim = max(pyclone_input_subset$mutation_count)
  
  
   # Rename legend for timepoints
  Timepoint <- factor(pyclone_input_subset$tumor_descriptor)
  
  # Plot stacked barplot 
  print(ggplot(pyclone_input_subset, aes(x = sample_id, 
                                              y = mutation_count, 
                                              fill = Timepoint)) +  
               geom_col(position = position_stack(reverse = TRUE)) +
               geom_bar(stat = "identity", width = 0.5) + 
               scale_fill_manual(values = palette, breaks = sort(names(palette))) + 
               theme_Publication() + 
               theme(axis.text.x = element_text(angle = 85, 
                                                hjust = 1, 
                                                vjust = 1)) + 
               labs(title = paste(sample_name)) + 
               labs(x = "sample_id", y = "Total Mutations") +
               ylim(0, ylim)) 
  
  
}
}
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KTRJ8TFY"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_KZ56XHJT"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"
[1] "PT_Z4BF2NSB"

sessionInfo()
R version 4.3.1 (2023-06-16)
Platform: x86_64-apple-darwin20 (64-bit)
Running under: macOS Ventura 13.6.1

Matrix products: default
BLAS:   /System/Library/Frameworks/Accelerate.framework/Versions/A/Frameworks/vecLib.framework/Versions/A/libBLAS.dylib 
LAPACK: /Library/Frameworks/R.framework/Versions/4.3-x86_64/Resources/lib/libRlapack.dylib;  LAPACK version 3.11.0

locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8

time zone: America/New_York
tzcode source: internal

attached base packages:
[1] grid      stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
 [1] ggthemes_4.2.4  readxl_1.4.3    lubridate_1.9.3 forcats_1.0.0   stringr_1.5.1   dplyr_1.1.4     purrr_1.0.2     readr_2.1.4     tidyr_1.3.0    
[10] tibble_3.2.1    ggplot2_3.4.4   tidyverse_2.0.0

loaded via a namespace (and not attached):
 [1] utf8_1.2.4        generics_0.1.3    stringi_1.8.1     hms_1.1.3         digest_0.6.33     magrittr_2.0.3    evaluate_0.23     timechange_0.2.0 
 [9] fastmap_1.1.1     cellranger_1.1.0  rprojroot_2.0.4   fansi_1.0.5       scales_1.2.1      clonevol_0.99.11  textshaping_0.3.7 cli_3.6.1        
[17] rlang_1.1.2       crayon_1.5.2      bit64_4.0.5       munsell_0.5.0     withr_2.5.2       yaml_2.3.7        tools_4.3.1       parallel_4.3.1   
[25] tzdb_0.4.0        colorspace_2.1-0  vctrs_0.6.4       R6_2.5.1          lifecycle_1.0.4   bit_4.0.5         vroom_1.6.4       ragg_1.2.6       
[33] pkgconfig_2.0.3   pillar_1.9.0      gtable_0.3.4      glue_1.6.2        systemfonts_1.0.5 xfun_0.41         tidyselect_1.2.0  rstudioapi_0.15.0
[41] knitr_1.45        farver_2.1.1      htmltools_0.5.7   rmarkdown_2.25    labeling_0.4.3    compiler_4.3.1   
LS0tCnRpdGxlOiAiSW5mZXJlbmNlIG9mIHN1YmNsb25hbCBhcmNoaXRlY3R1cmUgb2YgdHVtb3JzIGFjcm9zcyBtdWx0aXBsZSB0aW1lcG9pbnRzIGluIHRoZSBwYWlyZWQgbG9uZ2l0dWRpbmFsIChQTCkgY29ob3J0IgphdXRob3I6ICdBbnRvbmlhIENocm9uaSA8Y2hyb25pYUBjaG9wLmVkdT4gZm9yIEQzQicKZGF0ZTogIjIwMjMiCm91dHB1dDoKICBodG1sX25vdGVib29rOgogICAgdG9jOiBUUlVFCiAgICB0b2NfZmxvYXQ6IFRSVUUKLS0tCgojIEJhY2tncm91bmQKClRoaXMgbm90ZWJvb2sgaXMgZm9yIHByZXBhcmluZyBpbnB1dCBmaWxlcyBmb3IgW3B5Y2xvbmUtdmldKGh0dHBzOi8vZ2l0aHViLmNvbS9Sb3RoLUxhYi9weWNsb25lLXZpKS4gVGhlIHNhbWUgZmlsZXMgd2l0aCBhIG1pbm9yIG1vZGlmaWNhdGlvbiB3aWxsIGJlIHVzZWQgYXMgYW4gaW5wdXQgZm9yIFtGYXN0Q2xvbmVdKGh0dHBzOi8vZ2l0aHViLmNvbS9HdWFuTGFiL0Zhc3RDbG9uZV9HdWFuTGFiKS4KCgojIFNldCB1cAoKYGBge3IgbG9hZC1saWJyYXJ5fQpzdXBwcmVzc1BhY2thZ2VTdGFydHVwTWVzc2FnZXMoewogIGxpYnJhcnkodGlkeXZlcnNlKQogIGxpYnJhcnkocmVhZHhsKQogICNsaWJyYXJ5KGNEcml2ZXIpICMgQ2FsY3VsYXRlIENDRiwgaHR0cHM6Ly9naXRodWIuY29tL2hhbmFzdXNhay9jRHJpdmVyLwp9KQpgYGAKCiMgRGlyZWN0b3JpZXMgYW5kIEZpbGUgSW5wdXRzL091dHB1dHMKCmBgYHtyIHNldC1kaXItYW5kLWZpbGUtbmFtZXN9CiMgRGV0ZWN0IHRoZSAiLmdpdCIgZm9sZGVyIC0tIHRoaXMgd2lsbCBiZSBpbiB0aGUgcHJvamVjdCByb290IGRpcmVjdG9yeQojIFVzZSB0aGlzIGFzIHRoZSByb290IGRpcmVjdG9yeSB0byBlbnN1cmUgcHJvcGVyIHNvdXJjaW5nIG9mIGZ1bmN0aW9ucyAKIyBubyBtYXR0ZXIgd2hlcmUgdGhpcyBpcyBjYWxsZWQgZnJvbQpyb290X2RpciA8LSBycHJvanJvb3Q6OmZpbmRfcm9vdChycHJvanJvb3Q6Omhhc19kaXIoIi5naXQiKSkKYW5hbHlzaXNfZGlyIDwtIGZpbGUucGF0aChyb290X2RpciwgImFuYWx5c2VzIiwgInR1bW9yLWNsb25lLWluZmVyZW5jZSIpCmlucHV0X2RpciA8LSBmaWxlLnBhdGgoYW5hbHlzaXNfZGlyLCAiaW5wdXQiKQpkYXRhX2RpciA8LSBmaWxlLnBhdGgocm9vdF9kaXIsICJkYXRhIikKZmlsZXNfZGlyIDwtIGZpbGUucGF0aChyb290X2RpciwgImFuYWx5c2VzIiwgInNhbXBsZS1kaXN0cmlidXRpb24tYW5hbHlzaXMiLCAicmVzdWx0cyIpCm1hZl9maWxlc19kaXIgPC0gZmlsZS5wYXRoKHJvb3RfZGlyLCAiYW5hbHlzZXMiLCAidG1iLXZhZi1sb25naXR1ZGluYWwiLCAicmVzdWx0cyIpCiN0bWJfaW5wdXRfZGlyIDwtIGZpbGUucGF0aChyb290X2RpciwgImFuYWx5c2VzIiwgInRtYi12YWYtbG9uZ2l0dWRpbmFsIiwgImlucHV0IikKI3NjcmF0Y2hfZGlyIDwtIGZpbGUucGF0aChyb290X2RpciwgInNjcmF0Y2giKQoKIyBJbnB1dCBmaWxlcwpwYnRhX2ZpbGUgPC0gZmlsZS5wYXRoKGZpbGVzX2RpciwgInBidGEudHN2IikgIyBmaWxlIGZyb20gYWRkLXNhbXBsZS1kaXN0cmlidXRpb24gbW9kdWxlCm1hZl9maWxlIDwtIGZpbGUucGF0aChkYXRhX2RpciwgInNudi1jb25zZW5zdXMtcGx1cy1ob3RzcG90cy5tYWYudHN2Lmd6IikKdG1iX2ZpbGUgPC0gZmlsZS5wYXRoKG1hZl9maWxlc19kaXIsICJ0bWJfdmFmX2dlbm9taWMudHN2IikKZ2Vub21pY19wYWlyZWRfZmlsZSA8LSBmaWxlLnBhdGgoZmlsZXNfZGlyLCAiZ2Vub21pY19hc3NheXNfbWF0Y2hlZF90aW1lX3BvaW50cy50c3YiKSAjIGZpbGUgZnJvbSBhZGQtc2FtcGxlLWRpc3RyaWJ1dGlvbiBtb2R1bGUKbmF1dGlsdXNfZGVjX2ZpbGUgPC0gZmlsZS5wYXRoKGlucHV0X2RpciwgImRlY2Vhc2VkX3NhbXBsZXMueGxzeCIpIApwYWxldHRlX2ZpbGUgPC0gZmlsZS5wYXRoKHJvb3RfZGlyLCAiZmlndXJlcyIsICJwYWxldHRlcyIsICJ0dW1vcl9kZXNjcmlwdG9yX2NvbG9yX3BhbGV0dGUudHN2IikKI3RtYl9hbGxfZmlsZSA8LSBmaWxlLnBhdGgodG1iX2lucHV0X2RpciwgInNudi1tdXRhdGlvbi10bWItY29kaW5nLnRzdiIpCgoKIyBFWEFNUExFIAojIEJ1dCB0byByZXBsYWNlIHdpdGggZGlyIHdpdGggY25zIGRhdGEgaW5mZXJyZWQgYnkgQ05Wa2l0CiMgUFRfWjRCRjJOU0JfZGlyIDwtIGZpbGUucGF0aChpbnB1dF9kaXIsICJjbnZraXRfZGF0YV9leGFtcGxlL1BUX1o0QkYyTlNCIikgCmNuc19kaXIgPC0gZmlsZS5wYXRoKGlucHV0X2RpciwgImNudmtpdF9kYXRhIikgCgojIEZpbGUgcGF0aCB0byByZXN1bHRzIGRpcmVjdG9yeQpyZXN1bHRzX2RpciA8LQogIGZpbGUucGF0aChhbmFseXNpc19kaXIsICJyZXN1bHRzIikKaWYgKCFkaXIuZXhpc3RzKHJlc3VsdHNfZGlyKSkgewogIGRpci5jcmVhdGUocmVzdWx0c19kaXIpCn0KCiMgRmlsZSBwYXRoIHRvIGlucHV0IGRpcmVjdG9yeQpweWNsb25ldmlfaW5wdXRfZGlyIDwtCiAgZmlsZS5wYXRoKGFuYWx5c2lzX2RpciwgInJlc3VsdHMiLCAicHljbG9uZS12aS1pbnB1dCIpCmlmICghZGlyLmV4aXN0cyhweWNsb25ldmlfaW5wdXRfZGlyKSkgewogIGRpci5jcmVhdGUocHljbG9uZXZpX2lucHV0X2RpcikKfQoKIyBGaWxlIHBhdGggdG8gaW5wdXQgZGlyZWN0b3J5CmZhc3RjbG9uZV9pbnB1dF9kaXIgPC0KICBmaWxlLnBhdGgoYW5hbHlzaXNfZGlyLCAicmVzdWx0cyIsICJmYXN0Y2xvbmUtaW5wdXQiKQppZiAoIWRpci5leGlzdHMocHljbG9uZXZpX2lucHV0X2RpcikpIHsKICBkaXIuY3JlYXRlKHB5Y2xvbmV2aV9pbnB1dF9kaXIpCn0KCgojIEZpbGUgcGF0aCB0byBwbG90IGRpcmVjdG9yeQpweWNsb25lX3Bsb3RzX2RpciA8LQogIGZpbGUucGF0aChhbmFseXNpc19kaXIsICJwbG90cyIsICJweWNsb25lLXZpIikKaWYgKCFkaXIuZXhpc3RzKHB5Y2xvbmVfcGxvdHNfZGlyKSkgewogIGRpci5jcmVhdGUocHljbG9uZV9wbG90c19kaXIpCn0KCgpzb3VyY2UocGFzdGUwKHJvb3RfZGlyLCAiL2ZpZ3VyZXMvc2NyaXB0cy90aGVtZS5SIikpCmBgYAoKIyBMb2FkIGFuZCBwcm9jZXNzIGRhdGEKCmBgYHtyIGxvYWQtcHJvY2Vzcy1pbnB1dHN9CiMgUmVhZCBpbiBoaXN0b2xvZ2llcyBmaWxlIGFuZCBmaWx0ZXIgZm9yIHRoZSBwYnRhIGNvaG9ydApwYnRhIDwtIHJlYWRyOjpyZWFkX3RzdihwYnRhX2ZpbGUsIGd1ZXNzX21heCA9IDEwMDAwMCwgc2hvd19jb2xfdHlwZXMgPSBGQUxTRSkgJT4lCiAgZmlsdGVyKCEoZXhwZXJpbWVudGFsX3N0cmF0ZWd5ID09ICJSTkEtU2VxIikpICU+JSAKICBzZWxlY3QoS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCwgS2lkc19GaXJzdF9CaW9zcGVjaW1lbl9JRCwgY2dfaWQsIHR1bW9yX2Rlc2NyaXB0b3IsCiAgICAgICAgIHNhbXBsZV9pZCwgYWxpcXVvdF9pZCwKICAgICAgICAgdHVtb3JfZnJhY3Rpb24sIHR1bW9yX3Bsb2lkeSkKCiMgUmVhZCBpbiB0bWJfYWxsIGZpbGUKdG1iX2RmIDwtIHJlYWRyOjpyZWFkX3Rzdih0bWJfZmlsZSwgZ3Vlc3NfbWF4ID0gMTAwMDAwLCBzaG93X2NvbF90eXBlcyA9IEZBTFNFKSAKCiMgUmVhZCBtYWYKIyBLZWVwIGJvdGggc3lub255bW91cyBhbmQgbm9uLXN5bm95bW91cyBtdXRhdGlvbnMgZm9yIHRoZSBjdXJyZW50IGFuYWx5c2lzCm1hZl9maWxlX2RmIDwtIHZyb29tOjp2cm9vbShtYWZfZmlsZSwgZGVsaW0gPSAiXHQiLCBjb2xfbmFtZXMgPSBUUlVFKQoKIyBQcm9jZXNzIG1hZgptYWZfZGYgPC0gbWFmX2ZpbGVfZGYgJT4lCiAgZmlsdGVyKFR1bW9yX1NhbXBsZV9CYXJjb2RlICVpbiUgcGJ0YSRLaWRzX0ZpcnN0X0Jpb3NwZWNpbWVuX0lEIywKICAgICAgICAgI1ZhcmlhbnRfQ2xhc3NpZmljYXRpb24gJWluJSBtYWZfbm9uc3lub255bW91cwogICAgICAgICApICU+JQogIGRwbHlyOjpyZW5hbWUoS2lkc19GaXJzdF9CaW9zcGVjaW1lbl9JRCA9IFR1bW9yX1NhbXBsZV9CYXJjb2RlKSAlPiUgCiAgbGVmdF9qb2luKHBidGEpICU+JQogIGxlZnRfam9pbih0bWJfZGYpICU+JQoKICAjIENhbGN1bGF0ZSBWQUYgCiAgZHBseXI6Om11dGF0ZShWQUYgPSB0X2FsdF9jb3VudCAvICh0X3JlZl9jb3VudCArIHRfYWx0X2NvdW50KSwgCiAgICAgICAgICAgICAgICAjIENvbmNhdGVuYXRlIGdlbmUgYW5kIHByb3RlaW4gaW5mb3JtYXRpb24KICAgICAgICAgICAgICAgIGdlbmVfcHJvdGVpbiA9IHBhc3RlKEh1Z29fU3ltYm9sLCBIR1ZTcF9TaG9ydCwgc2VwID0gIl8iKSkgJT4lICAKICBzZWxlY3QoS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCwgS2lkc19GaXJzdF9CaW9zcGVjaW1lbl9JRCwgY2dfaWQsIHR1bW9yX2Rlc2NyaXB0b3IsCiAgICAgICAgIHNhbXBsZV9pZCwgYWxpcXVvdF9pZCwgQ2hyb21vc29tZSwgCiAgICAgICAgIFN0YXJ0X1Bvc2l0aW9uLCBSZWZlcmVuY2VfQWxsZWxlLCBUdW1vcl9TZXFfQWxsZWxlMSwgVHVtb3JfU2VxX0FsbGVsZTIsCiAgICAgICAgIHRfcmVmX2NvdW50LCB0X2FsdF9jb3VudCwgCiAgICAgICAgIHR1bW9yX2ZyYWN0aW9uLCB0dW1vcl9wbG9pZHksIFZBRiwgZ2VuZV9wcm90ZWluLCBIdWdvX1N5bWJvbCwKICAgICAgICAgbXV0YXRpb25fY291bnQsIHRtYikgJT4lIAoKICAjIFJlbW92ZSBoeXBlcm11dGFudHMKICBmaWx0ZXIoIXRtYiA+PSAxMCAjLCAKICAKICAjIFRoZXJlIGFyZSBhbHRlcmF0aW9ucyB3aXRoIGhpZ2ggbnVtYmVyIG9mIHJlYWQgY291bnRzLgogICMgV2Ugd2lsbCBleGNsdWRlIHRob3NlIHdpdGggPjEwMDAgZm9yIG5vdy4KICAgICAgICAjICF0X2FsdF9jb3VudCA+PSAxMDAwLAogICAgICAgICAjIXRfcmVmX2NvdW50ID49IDEwMDAKICApICAlPiUgCiAgCiAgIyBBZGQgYG5vcm1hbF9jbmA6IFRvdGFsIGNvcHkgbnVtYmVyIG9mIHNlZ21lbnQgaW4gaGVhbHRoeSB0aXNzdWUuIAogICMgRm9yIGF1dG9zb21lIHRoaXMgd2lsbCBiZSB0d28gYW5kIG1hbGUgc2V4IGNocm9tb3NvbWVzIG9uZS4KICAjIFNlZTogaHR0cHM6Ly9naXRodWIuY29tL1JvdGgtTGFiL3B5Y2xvbmUtdmkKICBtdXRhdGUobm9ybWFsX2NuID0gY2FzZV93aGVuKGdyZXBsKCJjaHJZIiwgQ2hyb21vc29tZSkgfiAiMSIsIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgVFJVRSB+ICIyIikpICU+JSAKICAKICAjIENhbGN1bGF0ZSBDQ0YKICAjIENDRigpICU+JSAKICAjIHVzaW5nIGZvcm11bGEgZm9yIENDRiBjYWxjdWxhdGlvbgogIG11dGF0ZShDQ0YgPSAoMip0X2FsdF9jb3VudCkvICh0X3JlZl9jb3VudCArIHRfYWx0X2NvdW50KSkgIyU+JSAKICAjd3JpdGVfdHN2KGZpbGUucGF0aChzY3JhdGNoX2RpciwgIm1hZl9wYnRhLnRzdiIpKQoKIyBOdW1iZXIgb2Ygc2FtcGxlcyBpbiB0aGUgYG1hZl9kZmAKbWFmX3NhbXBsZXMgPC0gcHJpbnQobGVuZ3RoKHVuaXF1ZShtYWZfZGYkS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCkpKQoKIyBSZWFkIHBhdGllbnQgbGlzdApnZW5vbWljX3BhaXJlZF9kZiA8LSByZWFkcjo6cmVhZF90c3YoZ2Vub21pY19wYWlyZWRfZmlsZSwgZ3Vlc3NfbWF4ID0gMTAwMDAwLCBzaG93X2NvbF90eXBlcyA9IEZBTFNFKSAlPiUgCiAgc2VsZWN0KCFjKGNhbmNlcl9ncm91cCwgZXhwZXJpbWVudGFsX3N0cmF0ZWd5KSkgJT4lIAogIGxlZnRfam9pbihtYWZfZGYsIGJ5ID0gYygiS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCIpKSAlPiUgIyBjcmVhdGUgdW5pcXVlIGlkZW50aWZpZXJzCiAgZHBseXI6Om11dGF0ZShtdXRhdGlvbl9pZCA9IHBhc3RlKEtpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQsIHR1bW9yX2Rlc2NyaXB0b3IsIEtpZHNfRmlyc3RfQmlvc3BlY2ltZW5fSUQsIENocm9tb3NvbWUsIFN0YXJ0X1Bvc2l0aW9uLCBSZWZlcmVuY2VfQWxsZWxlLCBUdW1vcl9TZXFfQWxsZWxlMiwgc2VwID0gIjoiKSwKICAgICAgICAgICAgICAgIGNnX2lkX2tpZHMgPSBwYXN0ZShjZ19pZCwgS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCwgc2VwID0gIl8iKSwKICAgICAgICAgICAgICAgIGNnX2lkX2tpZHMgPSBzdHJfcmVwbGFjZShjZ19pZF9raWRzLCAiL3wtIiwgIl8iKSwKICAgICAgICAgICAgICAgIGNnX2lkX2tpZHMgPSBzdHJfcmVwbGFjZV9hbGwoY2dfaWRfa2lkcywgIiAiLCAiXyIpLAogICAgICAgICAgICAgICAgY2dfaWQgPSBzdHJfcmVwbGFjZShjZ19pZCwgIi98LSIsICJfIiksCiAgICAgICAgICAgICAgICBjZ19pZCA9IHN0cl9yZXBsYWNlX2FsbChjZ19pZCwgIiAiLCAiXyIpLAogICAgICAgICAgICAgICAgc2FtcGxlX2lkID0gcGFzdGUodHVtb3JfZGVzY3JpcHRvciwgS2lkc19GaXJzdF9CaW9zcGVjaW1lbl9JRCwgc2VwID0gIjoiKSkgIyU+JSAKICAjZmlsdGVyKCFpcy5uYShWQUYpICMsCiAgICAgICAgICMhVkFGID09IDAgIywgCiAgICAgICAgICMhaXMubmEodG1iKQogICAgICAgICMgKQoKIyBOdW1iZXIgb2Ygc2FtcGxlcyBpbiB0aGUgYGdlbm9taWNfcGFpcmVkX2RmYApnZW5vbWljX3NhbXBsZXMgPC0gcHJpbnQobGVuZ3RoKHVuaXF1ZShnZW5vbWljX3BhaXJlZF9kZiRLaWRzX0ZpcnN0X1BhcnRpY2lwYW50X0lEKSkpCgojIExldCdzIGNvdW50ICNzcGVjaW1lbnMgcGVyIHNhbXBsZSAKIyBXZSB3aWxsIHJlbW92ZSBhbnkgc2FtcGxlcyB3aXRoIGxlc3MgdGhhbiAyIHNwZWNpbWVucyBhcyBwaHlsb2dlbmllcyByZXF1aXJlIGF0IGxlYXN0IDMgdGF4YQpraWRzX3NwZWNpbWVuc19uX2RmIDwtIGdlbm9taWNfcGFpcmVkX2RmICU+JSAKICBzZWxlY3QoS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCwgS2lkc19GaXJzdF9CaW9zcGVjaW1lbl9JRCwgdHVtb3JfZGVzY3JpcHRvcikgJT4lIAogIHVuaXF1ZSgpICU+JSAKICBkcGx5cjo6Y291bnQoS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCkgJT4lIAogIGRwbHlyOjptdXRhdGUoa2lkc19zcGVjaW1lbnNfbiA9IGdsdWU6OmdsdWUoIntLaWRzX0ZpcnN0X1BhcnRpY2lwYW50X0lEfSAgKE49e259KSIpKSAlPiUKICBkcGx5cjo6cmVuYW1lKGtpZHNfc3BlY2ltZW5zX251bWJlciA9IG4pICU+JSAKICBmaWx0ZXIoIWtpZHNfc3BlY2ltZW5zX251bWJlciA8PSAyKSAlPiUKCiAgbGVmdF9qb2luKGdlbm9taWNfcGFpcmVkX2RmLCBieSA9IGMoIktpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQiKSkgJT4lCiAgI2xlZnRfam9pbihtYWZfZGYpICU+JSAKICBmaWx0ZXIoIWlzLm5hKENocm9tb3NvbWUpKSAKCgojIExldCdzIGNvbmZpcm0gYW5kIGFkZCB0aGUgbnVtYmVyIG9mIHRpbWVwb2ludHMgcGVyIHNhbXBsZQojIEluIGNhc2Ugd2UgbG9zdCBhbnkgZnJvbSBmaWx0ZXJpbmcgaHlwZXJtdXRhbnRzIGFuZCBoaWdoIHJlYWRzL2FsdGVyYXRpb24KdGltZXBvaW50c19uX2RmIDwtIGtpZHNfc3BlY2ltZW5zX25fZGYgJT4lIAogIHNlbGVjdChLaWRzX0ZpcnN0X1BhcnRpY2lwYW50X0lELCB0dW1vcl9kZXNjcmlwdG9yKSAlPiUgCiAgdW5pcXVlKCkgJT4lIAogIGRwbHlyOjpjb3VudChLaWRzX0ZpcnN0X1BhcnRpY2lwYW50X0lEKSAlPiUgCiAgZHBseXI6Om11dGF0ZShraWRzX3RpbWVwb2ludHNfbiA9IGdsdWU6OmdsdWUoIntLaWRzX0ZpcnN0X1BhcnRpY2lwYW50X0lEfSAgKE49e259KSIpKSAlPiUKICBkcGx5cjo6cmVuYW1lKGtpZHNfdGltZXBvaW50c19udW1iZXIgPSBuKSAlPiUgCiAgZmlsdGVyKCFraWRzX3RpbWVwb2ludHNfbnVtYmVyIDw9IDEpCgpkZiA8LSB0aW1lcG9pbnRzX25fZGYgJT4lIAogIGxlZnRfam9pbihraWRzX3NwZWNpbWVuc19uX2RmLCBieSA9IGMoIktpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQiKSkgJT4lIAogIHdyaXRlX3RzdihmaWxlLnBhdGgocmVzdWx0c19kaXIsICJzYW1wbGVzX2VsaWdpYmxlX2Zvcl9waHlsb2dlbnkudHN2IikpCgojIE51bWJlciBvZiBzYW1wbGVzIGluIHRoZSBgZ2Vub21pY19wYWlyZWRfZGZgIHdpdGggPiAzIGJpb3NwZWNpbWVucyBhbmQgPjIgdGltZXBvaW50cwpzYW1wbGVzX3dpdGhfM2JzIDwtIHByaW50KGxlbmd0aCh1bmlxdWUoZGYkS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCkpKQpzYW1wbGVzX3dpdGhfM2JzIDwtIHByaW50KHVuaXF1ZShkZiRLaWRzX0ZpcnN0X1BhcnRpY2lwYW50X0lEKSkKICAKIyBMaXN0IHdpdGggc2FtcGxlcyBlbGlnaWJsZSBmb3IgcGh5bG9nZW55CiMgSSBhZGRlZCB0aGUgaW5mb3JtYXRpb24gYWJvdXQgdG8gdXNlIGZvciBwaHlsb2dlbmV0aWMgaW5mZXJlbmNlcwpsaXN0X2RmIDwtIGRmICU+JSAKICBzZWxlY3QoS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCkgJT4lIAogIHVuaXF1ZSgpICU+JSAKICBtdXRhdGUoc29tYXRpY19nZXJtbGluZV9waHlsb2dlbnkgPSBjYXNlX3doZW4oZ3JlcGwoIlBUX0taNTZYSEpUfFBUX0tUUko4VEZZIiwgS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCkgfiAieWVzIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIFRSVUUgfiAibm90X3lldCIpKSAlPiUgCiAgd3JpdGVfdHN2KGZpbGUucGF0aChyZXN1bHRzX2RpciwgInNhbXBsZXNfZWxpZ2libGVfZm9yX3BoeWxvZ2VueV9saXN0LnRzdiIpKQoKIyBSZW1vdmUgbGFyZ2UgZmlsZXMKcm0oZ2Vub21pY19wYWlyZWRfZGYsIGtpZHNfc3BlY2ltZW5zX25fZGYsIHRtYl9kZikKCmBgYAoKIyMgQWRkIE5hdXRpbHVzIGxvY2F0aW9uIGZvciBEZWNlYXNlZCBzcGVjaW1lbnMKCmBgYHtyIGFkZC1OYXV0aWx1cy1sb2NhdGlvbn0KbmF1dGlsdXNfZGVjX2RmIDwtIHJlYWRfZXhjZWwobmF1dGlsdXNfZGVjX2ZpbGUpICU+JSAKICByaWdodF9qb2luKGRmLCBieSA9IGMoInNhbXBsZV9pZCIsICJhbGlxdW90X2lkIiwgInR1bW9yX2Rlc2NyaXB0b3IiKSkgJT4lCiAgZHBseXI6Om11dGF0ZShtdXRhdGlvbl9pZCA9IHBhc3RlKEtpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQsIHR1bW9yX2Rlc2NyaXB0b3IsIEtpZHNfRmlyc3RfQmlvc3BlY2ltZW5fSUQsIENocm9tb3NvbWUsIFN0YXJ0X1Bvc2l0aW9uLCBSZWZlcmVuY2VfQWxsZWxlLCBUdW1vcl9TZXFfQWxsZWxlMiwgc2VwID0gIjoiKSkgJT4lIAogIHdyaXRlX3RzdihmaWxlLnBhdGgocmVzdWx0c19kaXIsICJuYXV0aWx1c19kZWMudHN2IikpCgpsaXN0X2RmIDwtIG5hdXRpbHVzX2RlY19kZiAlPiUgCiAgc2VsZWN0KEtpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQsIEtpZHNfRmlyc3RfQmlvc3BlY2ltZW5fSUQsIGBOb3RlIGZpZWxkIGZyb20gTmF1dGlsdXMgb2YgaW5pdGlhbCBwYXJlbnRgKSAlPiUgCiAgdW5pcXVlKCkgJT4lIAogIHdyaXRlX3RzdihmaWxlLnBhdGgocmVzdWx0c19kaXIsICJuYXV0aWx1c19kZWNfbGlzdC50c3YiKSkKCmBgYAoKIyMgT3RoZXIsIG1heWJlIHRvIGRlbGV0ZQoKYGBgIHtyIG90aGVyLXRoaW5nc30KCiMgTWFrZSBsaXN0IHdpdGggc2FtcGxlcyB3aXRoID4yIGJpb3NwZWNpbWVucyBhdCBEZWNlYXNlZCB0aW1lcG9pbnQKZGVjX25fZGYgPC0gZGYgJT4lIAogIGZpbHRlcih0dW1vcl9kZXNjcmlwdG9yID09ICJEZWNlYXNlZCIpICU+JSAKICBzZWxlY3QoS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCwgdHVtb3JfZGVzY3JpcHRvciwgS2lkc19GaXJzdF9CaW9zcGVjaW1lbl9JRCkgJT4lIAogIHVuaXF1ZSgpICU+JSAKICBkcGx5cjo6Y291bnQoS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCkgJT4lIAogIyBkcGx5cjo6bXV0YXRlKGtpZHNfZGVjX24gPSBnbHVlOjpnbHVlKCJ7S2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRH0gIChOPXtufSkiKSkgJT4lCiAgZHBseXI6OnJlbmFtZShraWRzX2RlY2Vhc2VkX2JzX251bWJlciA9IG4pICU+JSAKICBmaWx0ZXIoIWtpZHNfZGVjZWFzZWRfYnNfbnVtYmVyIDw9IDEpICU+JSAKICB3cml0ZV90c3YoZmlsZS5wYXRoKHJlc3VsdHNfZGlyLCAia2lkc19kZWNfbXVsdGlwbGVfYnNfbGlzdC50c3YiKSkKCiMjIy0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCiMgbGV0J3MgbG9vayBpbnRvIFBUXzNLTTlXOFM4CiNQVF8zS005VzhTOF9kZiA8LSBuYXV0aWx1c19kZWNfZGYgJT4lIAojICBmaWx0ZXIoS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCA9PSAiUFRfM0tNOVc4UzgiLAogIyAgICAgICAgS2lkc19GaXJzdF9CaW9zcGVjaW1lbl9JRCA9PSAiQlNfMk5RWFk1MjgiKSAlPiUgCiMgIHNlbGVjdChLaWRzX0ZpcnN0X1BhcnRpY2lwYW50X0lELCBLaWRzX0ZpcnN0X0Jpb3NwZWNpbWVuX0lELCBgTm90ZSBmaWVsZCBmcm9tIE5hdXRpbHVzIG9mIGluaXRpYWwgcGFyZW50YCkgJT4lIAojICB1bmlxdWUoKSAlPiUgCiMgIHdyaXRlX3RzdihmaWxlLnBhdGgocmVzdWx0c19kaXIsICJuYXV0aWx1c19kZWNfbGlzdC50c3YiKSkKICAKI2JzX2lkIDwtIHByaW50KHVuaXF1ZShQVF8zS005VzhTOF9kZiRLaWRzX0ZpcnN0X0Jpb3NwZWNpbWVuX0lEKSkKYGBgCgoKIyBDcmVhdGUgcHljbG9uZSBpbnB1dCBmaWxlcyAKCldlIG5lZWQgdG8gZ2VuZXJhdGUgdGhlIGlucHV0IGZpbGVzIGFjY29yZGluZyB0byB0aGUgbWV0aG9kJ3MgdGVtcGxhdGUuIApQaHlsb2dlbmV0aWMgbWV0aG9kcyByZXF1aXJlIGF0IGxlYXN0IDIgc2FtcGxlcyBwZXIgdHVtb3Igc2l0ZSAobXVsdGlyZWdpb25hbCBzYW1wbGluZyBwZXIgYW5hdG9taWNhbCBzaXRlKS4KSGVyZSwgd2Ugd2lsbCBjb25zaWRlciBraWRzIHNhbXBsZXMgd2l0aCBtb3JlIHRoYW4gMiB0aW1lcG9pbnRzIHdpdGggb25lIG9yIG1vcmUgYmlvc3BlY2ltZW5zLiAKV2Ugd2lsbCBjb21wYXJlIGxhdGVyIGRpZmZlcmVuY2VzIGluIHNhbXBsZXMgd2l0aCBzaW5nbGUgdnMgbXVsdGlwbGUgYmlvc3BlY2ltZW5zLgoKYGBge3IgY3JlYXRlLXB5Y2xvbmUtYWxsLXNhbXBsZXMtZGZ9CiMgQ3JlYXRlIHB5Y2xvbmUgZGYgZm9yIGFsbCBzYW1wbGVzCnB5Y2xvbmVfYWxsX3NhbXBsZXNfZGYgPC0gbmF1dGlsdXNfZGVjX2RmICU+JQogIHNlbGVjdChLaWRzX0ZpcnN0X1BhcnRpY2lwYW50X0lELCBLaWRzX0ZpcnN0X0Jpb3NwZWNpbWVuX0lELCBjZ19pZCwgCiAgICAgICAgIGNnX2lkX2tpZHMsIG11dGF0aW9uX2lkLCBzYW1wbGVfaWQsCiAgICAgICAgIHR1bW9yX2Rlc2NyaXB0b3IsIENocm9tb3NvbWUsIFN0YXJ0X1Bvc2l0aW9uLCAKICAgICAgICAgUmVmZXJlbmNlX0FsbGVsZSwgVHVtb3JfU2VxX0FsbGVsZTIsIHRfcmVmX2NvdW50LCB0X2FsdF9jb3VudCwKICAgICAgICAgbm9ybWFsX2NuLCB0dW1vcl9mcmFjdGlvbiwgbXV0YXRpb25fY291bnQpICU+JSAKICAKICAjIGNoYW5nZSBuYW1lcyB0byBtYXRjaCBpbnB1dCByZXF1aXJlbWVudHMKICBkcGx5cjo6cmVuYW1lKCJyZWZfY291bnRzIiA9ICJ0X3JlZl9jb3VudCIsIAogICAgICAgICAgICAgICAgImFsdF9jb3VudHMiID0gInRfYWx0X2NvdW50IiwKICAgICAgICAgICAgICAgICJ0dW1vdXJfY29udGVudCIgPSAidHVtb3JfZnJhY3Rpb24iKSAKCnNhbXBsZXNfcHljbG9uZSA8LSBwcmludCh1bmlxdWUocHljbG9uZV9hbGxfc2FtcGxlc19kZiRLaWRzX0ZpcnN0X1BhcnRpY2lwYW50X0lEKSkKCiMgUmVtb3ZlIGxhcmdlIGZpbGVzCiMgcm0obmF1dGlsdXNfZGVjX2RmKQoKIyBJIHdpbGwgdGVzdCBvbmUgSEdHIGRhdGFzZXQgZm9yIG5vdwojIFBUX1o0QkYyTlNCCgojUFRfWjRCRjJOU0JfREYgPC0gcHljbG9uZV9hbGxfc2FtcGxlc19kZiAlPiUgCiMgIGZpbHRlcihLaWRzX0ZpcnN0X1BhcnRpY2lwYW50X0lEID09ICJQVF9aNEJGMk5TQiIpCiAgICAgICAgICAgICAgICAgICAKYGBgCgojIyBBZGQgbWFqb3JfY24sIG1pbm9yX2NuIGNvbHVtbnMgZnJvbSBjbnMgZmlsZXMKCmBgYHtyIHByb2Nlc3MtY25zLWZpbGVzfQpkYXRhX2RpciA8LSBkaXIocGF0aCA9IGNuc19kaXIsICBwYXR0ZXJuID0gIi5jYWxsLmNucyIsIGZ1bGwubmFtZXMgPSBUUlVFLCByZWN1cnNpdmUgPSBUUlVFKQoKIyBDcmVhdGUgbGlzdCAKZGF0YV9saXN0IDwtIGxpc3QoKQoKZm9yIChpIGluIDE6bGVuZ3RoKGRhdGFfZGlyKSApIHsgCiAgCiAgIyBDcmVhdGUgc2FtcGxlX25hbWUKICBzYW1wbGVfbmFtZSA8LSAgdW5pcXVlKGFzLmNoYXJhY3Rlcihnc3ViKCIuY2FsbC5jbnMiLCAiIiwgc3RyX3NwbGl0X2ZpeGVkKGRhdGFfZGlyW2ldLCAiLyIsIDEzKVssMTFdKSkpCiAgI3NhbXBsZV9uYW1lIDwtICB1bmlxdWUoYXMuY2hhcmFjdGVyKGdzdWIoIi5jYWxsLmNucyIsICIiLCBzdHJfc3BsaXRfZml4ZWQoZGF0YV9kaXJbaV0sICIvIiwgMTApWyw5XSkpKQogIHNhbXBsZV9uYW1lIDwtIHNvcnQoc2FtcGxlX25hbWUsIGRlY3JlYXNpbmcgPSBGQUxTRSkKICBwcmludChzYW1wbGVfbmFtZSkKICAKICBmb3IgKHggaW4gc2VxX2Fsb25nKHNhbXBsZV9uYW1lKSApIHsgCiAgICAKICAgIGRhdGFfbGlzdFtbaV1dIDwtIHJlYWQuY3N2KGRhdGFfZGlyW2ldLCBoZWFkZXI9VCwgc2VwPSJcdCIpIAogIAogICAgCiAgIyBDcmVhdGUgZmlsZV9uYW1lCiAgZmlsZV9uYW1lIDwtIGdzdWIoIi5jYWxsLmNucyIsICIiLCBzdHJfc3BsaXRfZml4ZWQoZGF0YV9kaXJbaV0sICIvIiwgMTMpWywxMl0pCiAgcHJpbnQoZmlsZV9uYW1lKQogIAogIGRhdGFfbGlzdFtbaV1dIDwtIGRhdGFfbGlzdFtbaV1dICU+JSAKICAgIG11dGF0ZShLaWRzX0ZpcnN0X0Jpb3NwZWNpbWVuX0lEID0gZmlsZV9uYW1lKQogIAogICMgVGhlIGZvbGxvd2luZyBjb2RlIGFzc2lnbnMgbmFtZSB0byBkZgogICMgQnV0IGFkZHMgYW4gZXh0cmEgZGYgaW50byB0aGUgbGlzdCAtIGRvbnQga25vdyB5ZXQgd2h5IGJ1dCBsZXQncyByZW1vdmUsIGl0J3Mgbm90IG5lY2Vzc2FyeQogICNkZiA8LSBhc3NpZ24oZmlsZV9uYW1lLCBkYXRhX2xpc3QpIAogICNkYXRhX2xpc3RbW2ZpbGVfbmFtZV1dIDwtIGRmCiAgfQp9CgojIEJpbmQgYWxsIGRmIGZyb20gbGlzdCAKZGF0YV9saXN0X2JpbmQgPC0gZHBseXI6OmJpbmRfcm93cyhkYXRhX2xpc3QpIAoKIyBSZW1vdmUgbGFyZ2UgZmlsZXMKcm0oZGF0YV9saXN0KQoKIyBDcmVhdGUgYW5kIHNhdmUgcHljbG9uZV9pbnB1dCEKcHljbG9uZV9pbnB1dCA8LSBkYXRhX2xpc3RfYmluZCAlPiUgCiAgCiAgIyBSZW5hbWUgdG8gbWF0Y2ggaW5wdXQgZm9ybWF0CiAgIyBUbyBmaWd1cmUgb3V0IGlmIHRoZSBhc3NpZ25tZW50IGlzIGNvcnJlY3QKICBkcGx5cjo6cmVuYW1lKCJtYWpvcl9jbiIgPSAiY24xIiwgCiAgICAgICAgICAgICAgICAibWlub3JfY24iID0gImNuMiIpICU+JSAKCiAgbGVmdF9qb2luKHB5Y2xvbmVfYWxsX3NhbXBsZXNfZGYpICU+JSAKICAKICBmaWx0ZXIoIWlzLm5hKG1ham9yX2NuKSwKICAgICAgICAgIWlzLm5hKG1pbm9yX2NuKSkgJT4lCiAgI2RwbHlyOjptdXRhdGUobXV0YXRpb25faWQgPSBwYXN0ZShtdXRhdGlvbl9pZCwgbWFqb3JfY24sIG1pbm9yX2NuLCByZWZfY291bnRzLCBhbHRfY291bnRzLCBzZXAgPSAiOiIpKSAlPiUgCiAgZHBseXI6Om11dGF0ZShtdXRhdGlvbl9pZCA9IHBhc3RlKEtpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQsIHR1bW9yX2Rlc2NyaXB0b3IsIEtpZHNfRmlyc3RfQmlvc3BlY2ltZW5fSUQsIENocm9tb3NvbWUsIFN0YXJ0X1Bvc2l0aW9uLCBSZWZlcmVuY2VfQWxsZWxlLCBUdW1vcl9TZXFfQWxsZWxlMiwgc2VwID0gIjoiKSkgJT4lIAogIHNlbGVjdChLaWRzX0ZpcnN0X1BhcnRpY2lwYW50X0lELCB0dW1vcl9kZXNjcmlwdG9yLCBjZ19pZF9raWRzLCBtdXRhdGlvbl9pZCwgc2FtcGxlX2lkLCByZWZfY291bnRzLCAKICAgICAgICAgYWx0X2NvdW50cywgbm9ybWFsX2NuLCBtYWpvcl9jbiwgbWlub3JfY24sIHR1bW91cl9jb250ZW50LCBtdXRhdGlvbl9jb3VudCkgJT4lIAogIAogICMgVG8gZW5zdXJlIHRoZXJlIGFyZSBubyBkdXBsaWNhdGVkIGVudHJpZXMgaW4gdGhlIGRhdGFmcmFtZQogIGRpc3RpbmN0KCkgJT4lIAogIGZpbHRlcighaXMubmEoS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCkpIAoKIyBTYXZlIGRmIHdpdGggYWxsIGRhdGEgdG9nZXRoZXIKIyBXZSB3aWxsIG5lZWQgdGhpcyB0byBwcmVwYXJlIHRoZSBpbnB1dCBmaWVsIGZvciBDbG9uRXZvbApweWNsb25lX2lucHV0X2FsbCA8LSBuYXV0aWx1c19kZWNfZGYgJT4lCiAgc2VsZWN0KG11dGF0aW9uX2lkLCBWQUYsIENDRiwgZ2VuZV9wcm90ZWluKSAlPiUgCiAgcmlnaHRfam9pbihweWNsb25lX2lucHV0KSAlPiUgCiAgd3JpdGVfdHN2KGZpbGUucGF0aChyZXN1bHRzX2RpciwgInB5Y2xvbmVfaW5wdXRfYWxsX2RhdGEudHN2IikpCgoKc2FtcGxlc19weWNsb25lIDwtIHByaW50KHVuaXF1ZShweWNsb25lX2lucHV0JEtpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQpKQoKCiMgVG8gZmluZCB0aGUgcG9zaXRpb24gb2YgZHVwbGljYXRlIGVsZW1lbnRzIGluIHgsIHVzZSB0aGlzOgojIGR1cGxpY2F0ZWQocHljbG9uZV9pbnB1dCkKCiMgRXh0cmFjdCBkdXBsaWNhdGUgZWxlbWVudHM6CiMgcHljbG9uZV9pbnB1dFtkdXBsaWNhdGVkKHB5Y2xvbmVfaW5wdXQpXQoKYGBgCgojIFNhdmUgaW5wdXQgZmlsZXMKCmBgYHtyIHNhdmUtaW5wdXQtZmlsZXN9CmZvciAoaSBpbiAxOmxlbmd0aChkYXRhX2RpcikgKSB7IAogIAogICMgQ3JlYXRlIHNhbXBsZV9uYW1lCiAgc2FtcGxlX25hbWUgPC0gIHVuaXF1ZShhcy5jaGFyYWN0ZXIoZ3N1YigiLmNhbGwuY25zIiwgIiIsIHN0cl9zcGxpdF9maXhlZChkYXRhX2RpcltpXSwgIi8iLCAxMylbLDExXSkpKQogICNzYW1wbGVfbmFtZSA8LSAgdW5pcXVlKGFzLmNoYXJhY3Rlcihnc3ViKCIuY2FsbC5jbnMiLCAiIiwgc3RyX3NwbGl0X2ZpeGVkKGRhdGFfZGlyW2ldLCAiLyIsIDEwKVssOV0pKSkKICBzYW1wbGVfbmFtZSA8LSBzb3J0KHNhbXBsZV9uYW1lLCBkZWNyZWFzaW5nID0gRkFMU0UpCiAgcHJpbnQoc2FtcGxlX25hbWUpCiAgCmZvciAoeCBpbiBzZXFfYWxvbmcoc2FtcGxlX25hbWUpICkgeyAKICAKICAjIFNhdmUgaW5wdXQgZmlsZSBmb3IgcHljbG9uZS12aQogIHB5Y2xvbmVfZm5hbWUgPC0gcGFzdGUwKHB5Y2xvbmV2aV9pbnB1dF9kaXIsICIvIiwgc2FtcGxlX25hbWVbeF0sICIudHN2IikKICBwcmludChweWNsb25lX2ZuYW1lKQogIAogIHB5Y2xvbmVfaW5wdXRfc3Vic2V0IDwtIHB5Y2xvbmVfaW5wdXQgJT4lCiAgICBmaWx0ZXIoS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCA9PSBzYW1wbGVfbmFtZVt4XSkgJT4lIAogICAgI3NlbGVjdCgtYyhLaWRzX0ZpcnN0X1BhcnRpY2lwYW50X0lEKSkgJT4lIAogIHdyaXRlX3RzdihmaWxlLnBhdGgocHljbG9uZV9mbmFtZSkpCiAgCiAgCiAgIyBTYXZlIGlucHV0IGZpbGUgZm9yIEZhc3RDbG9uZQogIGZhc3RjbG9uZV9mbmFtZSA8LSBwYXN0ZTAoZmFzdGNsb25lX2lucHV0X2RpciwgIi8iLCBzYW1wbGVfbmFtZVt4XSwgIi50c3YiKQogIHByaW50KGZhc3RjbG9uZV9mbmFtZSkKICAKICBmYXN0Y2xvbmVfaW5wdXRfc3Vic2V0IDwtIHB5Y2xvbmVfaW5wdXQgJT4lCiAgICBkcGx5cjo6cmVuYW1lKCJ2YXJfY291bnRzIiA9ICJhbHRfY291bnRzIikgJT4lIAogICAgZmlsdGVyKEtpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQgPT0gc2FtcGxlX25hbWVbeF0pICU+JSAKICAgICNzZWxlY3QoLWMoS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCkpICU+JSAKICB3cml0ZV90c3YoZmlsZS5wYXRoKGZhc3RjbG9uZV9mbmFtZSkpCiAgCiAgfQp9CgpgYGAKCiMgUGxvdCBkZXB0aCBjb3ZlcmFnZQoKYGBge3IgZGVmaW5lLXBhcmFtZXRlcnMtZm9yLXBsb3RzfQojIFJlYWQgY29sb3IgcGFsZXR0ZQpwYWxldHRlX2RmIDwtIHJlYWRyOjpyZWFkX3RzdihwYWxldHRlX2ZpbGUsIGd1ZXNzX21heCA9IDEwMDAwMCwgc2hvd19jb2xfdHlwZXMgPSBGQUxTRSkgJT4lIAogIG11dGF0ZSh0dW1vcl9kZXNjcmlwdG9yID0gY29sb3JfbmFtZXMpCgojIERlZmluZSBhbmQgb3JkZXIgcGFsZXR0ZQpwYWxldHRlIDwtIHBhbGV0dGVfZGYkaGV4X2NvZGVzCm5hbWVzKHBhbGV0dGUpIDwtIHBhbGV0dGVfZGYkdHVtb3JfZGVzY3JpcHRvcgoKIyBEZWZpbmUgdGltZXBvaW50cwp0aW1lcG9pbnRzID0gYygiRGlhZ25vc2lzIiwgIlByb2dyZXNzaXZlIiwgIlJlY3VycmVuY2UiLCAiRGVjZWFzZWQiLCAiU2Vjb25kIE1hbGlnbmFuY3kiLCAiVW5hdmFpbGFibGUiKQoKYGBgCgpgYGAge3IgZGVwdGgtY292ZXJhZ2UsIGZpZy53aWR0aCA9IDYsIGZpZy5oZWlnaHQgPSA1LCBmaWcuZnVsbHdpZHRoID0gVFJVRX0KZm9yIChpIGluIDE6bGVuZ3RoKGRhdGFfZGlyKSApIHsgCiAgCiAgIyBDcmVhdGUgc2FtcGxlX25hbWUKICBzYW1wbGVfbmFtZSA8LSAgdW5pcXVlKGFzLmNoYXJhY3Rlcihnc3ViKCIuY2FsbC5jbnMiLCAiIiwgc3RyX3NwbGl0X2ZpeGVkKGRhdGFfZGlyW2ldLCAiLyIsIDEzKVssMTFdKSkpCiAgc2FtcGxlX25hbWUgPC0gc29ydChzYW1wbGVfbmFtZSwgZGVjcmVhc2luZyA9IEZBTFNFKQogIHByaW50KHNhbXBsZV9uYW1lKQogIApmb3IgKHggaW4gc2VxX2Fsb25nKHNhbXBsZV9uYW1lKSApIHsgCiAgCiAgcHljbG9uZV9pbnB1dF9zdWJzZXQgPC0gcHljbG9uZV9pbnB1dCAlPiUKICAgIGZpbHRlcihLaWRzX0ZpcnN0X1BhcnRpY2lwYW50X0lEID09IHNhbXBsZV9uYW1lW3hdKSAlPiUgCiAgICBzZWxlY3QoLWMoS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCkpICMlPiUgCiAgICMgbXV0YXRlKHR1bW9yX2Rlc2NyaXB0b3IgPSBmYWN0b3IodHVtb3JfZGVzY3JpcHRvciksCiAgICAjICAgICB0dW1vcl9kZXNjcmlwdG9yID0gZmN0X3JlbGV2ZWwodHVtb3JfZGVzY3JpcHRvciwgdGltZXBvaW50cykpICU+JSAKICAgICNhcnJhbmdlKHR1bW9yX2Rlc2NyaXB0b3IsIHNhbXBsZV9pZCkKICAKICAjIE1ha2UgdGhpcyByZXByb2R1Y2libGUKICBzZXQuc2VlZCgyMDIzKQoKICAjIERlZmluZSBsYWJlbCBmb3IgcGxvdHMKICBUaW1lcG9pbnQgPC0gZmFjdG9yKHggPSBweWNsb25lX2lucHV0X3N1YnNldCR0dW1vcl9kZXNjcmlwdG9yLCBsZXZlbHMgPSB0aW1lcG9pbnRzKQogIAogICMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjCiAgIyBDcmVhdGUgYnhwIHJlZl9jb3VudHMKICAgIHAgPC0gcHJpbnQoZ2dwbG90KHB5Y2xvbmVfaW5wdXRfc3Vic2V0LCBhZXMoc2FtcGxlX2lkLCByZWZfY291bnRzLCBjb2xvciA9IFRpbWVwb2ludCkpICsgCiAgICAgICAgICAgICAgICAgZ2VvbV9qaXR0ZXIod2lkdGggPSAwLjE1LCBzaXplID0gMC43LCBhbHBoYSA9IDAuNikgKwogICAgICAgICAgICAgICAgIGdncGxvdDI6Omdlb21fYm94cGxvdChjb2xvciA9ICJibGFjayIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHNpemUgPSAwLjI1LAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICBhbHBoYSA9IDAsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGNvZWYgPSAwKSArICMgcmVtb3ZlIHdoaXNrZXJzCiAgICAgICAgICAgICAgICAgdGhlbWVfUHVibGljYXRpb24oKSArIAogICAgICAgICAgICAgICAgIHNjYWxlX2NvbG9yX21hbnVhbCh2YWx1ZXMgPSBwYWxldHRlLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBicmVha3MgPSBzb3J0KG5hbWVzKHBhbGV0dGUpKSkgKwogICAgICAgICAgICAgICAgICNyb3RhdGUoKSArCiAgICAgICAgICAgICAgICAgdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA5MCkpICsKICAgICAgICAgICAgICAgICBzdGF0X3N1bW1hcnkoZnVuLnk9bWVhbixzaGFwZT0xLGNvbD0nYmxhY2snLGdlb209J3BvaW50JykgKwogICAgICAgICAgICAgICAgIGxhYnModGl0bGUgPSBzYW1wbGVfbmFtZVt4XSwKICAgICAgICAgICAgICAgICAgICAgIHggPSAic2FtcGxlX2lkIiwKICAgICAgICAgICAgICAgICAgICAgIHkgPSAicmVmX2NvdW50cyIsCiAgICAgICAgICAgICAgICAgICAgICBjb2xvciA9ICJUaW1lcG9pbnQiKSkKICAgIAogICAgIyBTYXZlIHRoZSBwbG90CiAgICBnZ3NhdmUoZmlsZW5hbWUgPSBwYXN0ZTAoc2FtcGxlX25hbWVbeF0sICItcmVmX2NvdW50cy5wZGYiKSwgCiAgICAgICAgIHBhdGggPSBweWNsb25lX3Bsb3RzX2RpciwgCiAgICAgICAgIHdpZHRoID0gNiwgCiAgICAgICAgIGhlaWdodCA9IDUsIAogICAgICAgICBkZXZpY2UgPSAicGRmIiwgCiAgICAgICAgdXNlRGluZ2JhdHMgPSBGQUxTRSkKICAgIAoKICAjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIwogICMgQ3JlYXRlIGJ4cCBhbHRfY291bnRzCiAgICBwIDwtIHByaW50KGdncGxvdChweWNsb25lX2lucHV0X3N1YnNldCwgYWVzKHNhbXBsZV9pZCwgYWx0X2NvdW50cywgY29sb3IgPSBUaW1lcG9pbnQpKSArIAogICAgICAgICAgICAgICAgIGdlb21faml0dGVyKHdpZHRoID0gMC4xNSwgc2l6ZSA9IDAuNywgYWxwaGEgPSAwLjYpICsKICAgICAgICAgICAgICAgICBnZ3Bsb3QyOjpnZW9tX2JveHBsb3QoY29sb3IgPSAiYmxhY2siLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICBzaXplID0gMC4yNSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgYWxwaGEgPSAwLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICBjb2VmID0gMCkgKyAjIHJlbW92ZSB3aGlza2VycwogICAgICAgICAgICAgICAgIHRoZW1lX1B1YmxpY2F0aW9uKCkgKyAKICAgICAgICAgICAgICAgICBzY2FsZV9jb2xvcl9tYW51YWwodmFsdWVzID0gcGFsZXR0ZSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgYnJlYWtzID0gc29ydChuYW1lcyhwYWxldHRlKSkpICsKICAgICAgICAgICAgICAgICAjcm90YXRlKCkgKwogICAgICAgICAgICAgICAgIHRoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gOTApKSArCiAgICAgICAgICAgICAgICAgc3RhdF9zdW1tYXJ5KGZ1bi55PW1lYW4sc2hhcGU9MSxjb2w9J2JsYWNrJyxnZW9tPSdwb2ludCcpICsKICAgICAgICAgICAgICAgICBsYWJzKHRpdGxlID0gc2FtcGxlX25hbWVbeF0sCiAgICAgICAgICAgICAgICAgICAgICB4ID0gInNhbXBsZV9pZCIsCiAgICAgICAgICAgICAgICAgICAgICB5ID0gImFsdF9jb3VudHMiLAogICAgICAgICAgICAgICAgICAgICAgY29sb3IgPSAiVGltZXBvaW50IikpCiAgICAKICAgICMgU2F2ZSB0aGUgcGxvdAogICAgZ2dzYXZlKGZpbGVuYW1lID0gcGFzdGUwKHNhbXBsZV9uYW1lW3hdLCAiLWFsdF9jb3VudHMucGRmIiksIAogICAgICAgICBwYXRoID0gcHljbG9uZV9wbG90c19kaXIsIAogICAgICAgICB3aWR0aCA9IDYsIAogICAgICAgICBoZWlnaHQgPSA1LCAKICAgICAgICAgZGV2aWNlID0gInBkZiIsIAogICAgICAgIHVzZURpbmdiYXRzID0gRkFMU0UpCiAgCiAgCn0KfQoKYGBgCgoKIyBUb3RhbCBudW1iZXIgb2YgbXV0YXRpb25zIGFjcm9zcyB0aW1lcG9pbnRzIGFuZCBiaW9zcGVjaW1lbiBzYW1wbGUgcGVyIFBhdGllbnQgY2FzZQoKYGBge3IgY3JlYXRlLWJhcnBsb3Qtc2FtcGxlLCBmaWcud2lkdGggPSA2LCBmaWcuaGVpZ2h0ID0gNSwgZmlnLmZ1bGx3aWR0aCA9IFRSVUV9CmZvciAoaSBpbiAxOmxlbmd0aChkYXRhX2RpcikgKSB7IAogIAogICMgQ3JlYXRlIHNhbXBsZV9uYW1lCiAgc2FtcGxlX25hbWUgPC0gIHVuaXF1ZShhcy5jaGFyYWN0ZXIoZ3N1YigiLmNhbGwuY25zIiwgIiIsIHN0cl9zcGxpdF9maXhlZChkYXRhX2RpcltpXSwgIi8iLCAxMylbLDExXSkpKQogIHNhbXBsZV9uYW1lIDwtIHNvcnQoc2FtcGxlX25hbWUsIGRlY3JlYXNpbmcgPSBGQUxTRSkKICBwcmludChzYW1wbGVfbmFtZSkKICAKZm9yICh4IGluIHNlcV9hbG9uZyhzYW1wbGVfbmFtZSkgKSB7IAogIAogIHB5Y2xvbmVfaW5wdXRfc3Vic2V0IDwtIHB5Y2xvbmVfaW5wdXQgJT4lCiAgICBmaWx0ZXIoS2lkc19GaXJzdF9QYXJ0aWNpcGFudF9JRCA9PSBzYW1wbGVfbmFtZVt4XSkgJT4lIAogICAgc2VsZWN0KC1jKEtpZHNfRmlyc3RfUGFydGljaXBhbnRfSUQpKQogIAogICMgRGVmaW5lIHBhcmFtZXRlcnMgZm9yIGZ1bmN0aW9uCiAgeWxpbSA9IG1heChweWNsb25lX2lucHV0X3N1YnNldCRtdXRhdGlvbl9jb3VudCkKICAKICAKICAgIyBSZW5hbWUgbGVnZW5kIGZvciB0aW1lcG9pbnRzCiAgVGltZXBvaW50IDwtIGZhY3RvcihweWNsb25lX2lucHV0X3N1YnNldCR0dW1vcl9kZXNjcmlwdG9yKQogIAogICMgUGxvdCBzdGFja2VkIGJhcnBsb3QgCiAgcHJpbnQoZ2dwbG90KHB5Y2xvbmVfaW5wdXRfc3Vic2V0LCBhZXMoeCA9IHNhbXBsZV9pZCwgCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICB5ID0gbXV0YXRpb25fY291bnQsIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgZmlsbCA9IFRpbWVwb2ludCkpICsgIAogICAgICAgICAgICAgICBnZW9tX2NvbChwb3NpdGlvbiA9IHBvc2l0aW9uX3N0YWNrKHJldmVyc2UgPSBUUlVFKSkgKwogICAgICAgICAgICAgICBnZW9tX2JhcihzdGF0ID0gImlkZW50aXR5Iiwgd2lkdGggPSAwLjUpICsgCiAgICAgICAgICAgICAgIHNjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcyA9IHBhbGV0dGUsIGJyZWFrcyA9IHNvcnQobmFtZXMocGFsZXR0ZSkpKSArIAogICAgICAgICAgICAgICB0aGVtZV9QdWJsaWNhdGlvbigpICsgCiAgICAgICAgICAgICAgIHRoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gODUsIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBoanVzdCA9IDEsIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICB2anVzdCA9IDEpKSArIAogICAgICAgICAgICAgICBsYWJzKHRpdGxlID0gcGFzdGUoc2FtcGxlX25hbWUpKSArIAogICAgICAgICAgICAgICBsYWJzKHggPSAic2FtcGxlX2lkIiwgeSA9ICJUb3RhbCBNdXRhdGlvbnMiKSArCiAgICAgICAgICAgICAgIHlsaW0oMCwgeWxpbSkpIAogIAogIAp9Cn0KYGBgCgpgYGB7cn0Kc2Vzc2lvbkluZm8oKQpgYGAKCgo=